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1. INTRODUCTION 

Accurate measurements of oscillation frequencies of 
the Sun are providing detailed information about 
the structure of the Sun. These frequencies can 
be used to test the basic input physics, like opaci- 
ties, equation of state and nuclear energy generation 
rates which are required to construct theoretical so- 
lar models. Here we concentrate on the effects of 
opacity. Throughout we consider models calibrated 
to solar radius and luminosity, by adjusting the ini- 
tial composition and the mixing-length parameter. 

2. METHOD 
We consider opacity modifications given by 

log* = log k + /(T,T ) , (1) 

where kq is the unmodified opacity, obtained as a 
function of T, density p and composition and the 
function /(T, T ) has the form 



f(T,T ) = Aexp 



log T — log T 



(2) 



where log is logarithm to base 10. The constants A 
and A set the magnitude and width of the opacity 
modification. By varying To, we can investigate the 
effects of changes in different parts of the models. 

Assuming a linear relation between small modifi- 
cations in the opacity and the response of the model 
(see Tripathy & Christensen-Dalsgaard 1995), the ef- 
fect of an opacity change 5 log k(T) for any quantity 
F can be approximated by the relation 



SF f 

— = J K F (T)S log K(T)dlogT 



(3) 



The kernel Kp(T) may be estimated from the change 
5F corresponding to the opacity modification given 
in equation (1) for sufficiently small A and A. 

For this work we consider F = c(r) 2 , the squared 
sound speed; thus Sc 2 /c 2 is the relative squared 
sound-speed difference between the Sun or a test 
model and the reference model. This quantity can 
be obtained by a straightforward inversion of so- 
lar oscillation frequencies (e.g., Basu et ai, 1996). 
It can be related to an intrinsic opacity difference 
Slogn through equation (0), with a kernel K c (r,T) 
which may be estimated by applying modifications 
of the form given in equation (2) for a range of To. 
If the sound-speed difference between the Sun and 



the model is assumed to arise solely from opacity er- 
rors, these may then be determined by inversion of 
the relation between SlogK and Sc 2 /c 2 . 

Here we expand the opacity correction in terms of 
the basis functions given by equation (2), as 



JY 



Slog K (T)=J2bif(T,T i ) 



(4) 



The resulting change in the squared sound speed can 
be written as 



N 



' 1=1 



(5) 



where ipi(r) is the change (Sc 2 /c 2 )(r) resulting from 
applying the opacity change SlogK — /(T, T). It 
was determined by carrying out a solar evolution 
calculation with the modified opacity. Having thus 
obtained ipi(r) for a suitable set of T, equation (||) 
may be fitted to the actual Sc 2 /c 2 by means of a 
regularized least-squares fit to determine the coef- 
ficients bi, and hence the opacity correction SlogK. 
We have calculated ijji{r) at 57 values of T, with 
logT = 6.2,. ..,7.19, A = 0.02 and A = 0.02. 

The procedure was tested by applying it to differ- 
ences between two models: the reference model was 
computed with the Cox & Tabor (1976) opacities, 
while the test model used the Los Alamos Opacity 
Library (Huebner et ai. 1977). Frequency differences 
between the two models, for the mode set used in 
the analysis of solar data, were inverted to infer the 
sound-speed differences, and the result was fitted to 
obtain the corresponding opacity change. The result- 
ing inferred S log k was in good agreement with the 
actual intrinsic opacity difference between the two 
tables (at fixed T, density p and composition). Thus 
we conclude that the method allows a determination 
of the opacity correction required to match a given 
sound-speed difference. 

3. RESULTS FOR SOLAR DATA 

We have applied the procedure to sound-speed dif- 
ferences between the Sun and a solar model. The 
data, model and inversion results are described by 
Basu et ai, (1996). The data are a combination of 
low-degree BiSON data and LOWL data of degrees 
between 3 and 99. The reference model (Model S of 
Christensen-Dalsgaard et al. 1996) used OPAL opac- 
ities (e.g. Iglesias et al. 1992), and included settling 
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Figure 1: Comparison between the inferred difference 
in squared sound speed between the Sun and the model 
(Model S of Christensen-Dalsgaard et al. 1996), shown 
with symbols, and the model change resulting from an 
evolution calculation where the opacity change in Fig. 2 
was added to the opacity tables. The differences are 
shown in the sense (Sun) - (reference model) . 




Figure 2: Opacity difference fitted from inferred 
sound-speed difference between the Sun and Model S of 
Christensen-Dalsgaard et al. (1996). 

of helium and heavy elements. The inversion was 
carried out by means of the Subtractive Optimally 
Localized Averaging technique. The inferred differ- 
ences in squared sound speed are shown in Figure 
1. 

The opacity modification resulting from a fit to the 
sound-speed difference is shown in Fig. 2. This indi- 
cates that a sharp feature in Slogn is required just 
beneath the convection zone, with a second peak in 
the core. To test that this opacity modification does 
in fact reproduce the inferred sound-speed difference, 
we have subsequently recomputed the solar evolution 
model, adding 5 log k(T) to the opacity as interpo- 
lated from the tables. In Fig. 1, the solid line shows 
the resulting difference in squared sound speed, rel- 
ative to the reference model. It is evident that this 
is in fact quite close to the sound-speed differences 
inferred from the inversion. 



We have attempted to infer the corrections to the 
opacity required to match the sound speed in a state- 
of-the-art solar model to that inferred from hclioscis- 
mic inversion. We find that the sound-speed differ- 
ence can to a large extent be reproduced by a change 
in the opacity in the radiative interior. The required 
change, of order of a few percent, is probably within 
the general level of uncertainty of current opacity 
calculations. Since the energy transport in the so- 
lar convection zone is independent of opacity, the 
sound-speed difference does not give us any informa- 
tion about the opacity in that region. 

It must be noted that the sound-speed difference 
between two models is essentially insensitive to a 
change in the opacity by a constant factor within 
the radiative interior. The change in opacity is 
largely compensated by a change in the composition, 
the sound-speed profile remaining approximately the 
same. Thus in the fit in Fig. 2, <51ogK is determined 
only to within a constant. The freedom in the level 
of the relative opacity change could, for example, be 
utilized to choose a solution for which the change in 
the initial hydrogen abundance is small. 

We also note that rather similar effects on the 
sound speed can be obtained through suitable weak 
mixing of the solar interior or at the base of the con- 
vection zone (Basu & Antia, 1994). Therefore, it can- 
not be claimed that the opacity errors are the only, 
or even the dominant, actual error in current solar- 
model calculations. A more detailed investigation 
of the individual contributions to the opacity, and 
their likely uncertainty, is required to test whether 
the specific shape of 5 log k is realistic. 
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